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the double- logarithmic approximation of perturbative QCD in the region i < 1. Double 
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non-ladder diagrams. 
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1 Introduction 



The investigation of the structure functions g\ and #2 provides the basis for the theoretical 
description of polarization effects in deep inelastic lepton nucleon scattering. Particular in- 
terest has been given to the behavior of g% at small x: experimental data Jl| on both gf and 
g^ are limited to x > 10~ 2 , and numerical values of the moments T P,N = J Q dxg l ' therefore 
depend upon the extrapolation in the region x < 10~ 2 . 

In QCD the Q 2 evolution at not too small x of g\ is, like in the unpolarized case, de- 
scribed by the GLAP-evolution equations 0, [|. The prediction for the small-x behavior is 
that both the quark and the gluon polarized structure functions go as 

Aqfa Q 2 ), Agfa Q 2 ) ~ exp y 'const ■ a s ln(Q 2 //x 2 ) ln(l/x). (1.1) 

In particular, in the flavor singlet part the gluons and the quarks mix. In a recent paper 
0] it has been shown, for the flavor nonsinglet contribution to g\, that this simple double 
log extrapolation (1) of the GLAP evolution equations, in fact, strongly underestimates the 
rise at small x. The reason is that, at x -C 1, new double logarithmic contributions appear 
which are beyond the control of the standard evolution framework. To be more precise, in 
the region of not too small x where the standard analysis applies the leading behavior in the 
n-th order of a s is of the type 

(a, ln(Q 2 //i 2 ))> n (ln(l/a;)r + a n _ 1 (ln(l/x)^ 1 + ...] (1.2) 

(where \x denotes the renormalization scale, and for simplicity we restrict ourselves to a fixed 
a s ), whereas at very small x the dominant contributions are of the form: 

(a s ln(l/x)ln(l/x)) n . (1.3) 

These double logarithmic contributions are not included in the standard evolution scheme 
0. In H the sum of the double logarithms 

(a s ) n [b n (ln(l/x)) 2n + b n ^{Hl/x)) 2n - 1 ln(Q 2 //i 2 ) + ... + 6 (ln(l/x))"(ln(g 2 / /U 2 ))™ (1.4) 

has been shown to give rise to a power-like increase of the flavor nonsinglet structure function 
at small-x which is stronger than the GLAP predicition (1). 

In this paper we continue our investigation of the small-x behavior of g\ in the double 
logarithmic approximation, by calculating the flavour singlet contribution. Like in the fla- 
vor nonsinglet case, we derive and solve evolution equations which describe the dependence 
upon the infrared cutoff of the transverse momentum integrations. Our main results is the 
power-like growth of g\ at small-x (eqs.(4.21) - (4.23)). A discussion of the phenomenological 
implications of our result will be presented in a forthcoming paper. 
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2 Singlet contribution to g\ in lowest orders a s 

Let us consider the scattering process where the off-shell photon with momentum q (—q 2 = 
Q 2 > 0) strikes the polarized quark with momentum p. We assume that the total energy is 
much greater than the mass of the quark,m. 

s = 2pq » m 2 . (2.1) 

We follow the notation of ||, [| and introduce the scattering amplitudes Tj in the following 

way: 

T^ u = iJ d 4 xe iqx < N\T(J a (x)J u (0)) \N > 

= {-9uu + — 2-)7i + (Pu - q»— ){Pu - q v — )T 2 
qt q z 

+ie, ua pq a sP-T 3 + it^q a {s^pq) - p^sq))-T 4 . (2.2) 
pq pq 

The structure function g x is defined as 

9i = ~^-ImT 3 . (2.3) 



In the Born approximation the matrix element for this process is given by the Feynman 
graphs in Fig 1. Let us denote by Ml(Mr) the matrix element for the case when the initial 
quark has left (the right) helicity. Then we obtain 

H/fBorn ]\/fBorn „ „ 
T Born _ M L M R _ 2/ S s \ (? A] 

so that 

#f orn = -—ImM Born . (2.5) 
2n 

e q here stands for the electric charge of the initial quark. Expession ( |2.4| ) shows that T 3 Born 
is antisymmetrical with respect to replacing s by — s, i.e. in terms of the Regge theory 
rpBom j-^g negative/odd signature. QCD radiative corrections turn T 3 Born into T 3 . As usual, 
it consists of the flavor nonsinglet part T 3 NS and the flavour singlet part T 3 5 . The small-x 
behavior of the nonsinglet contribution has been discussed in M, and in this paper we 
shall consider the singlet part. All Feynman graphs contributing to the double-logarithmic 
approximation (DLA) of the nonsinglet part T 3 NS can be obtained from the graphs in Fig.l by 
adding the gluon lines (further rungs and nonladder bremsstrahlungs gluons). An example 
is shown in Fig. 2. In the present case of the flavor singlet we will have to include more 
general contributions, in particular gluon ladders. An example is shown in Fig. 3. As a 
result, we first have to consider the four different types of kernels: AP qq , AP qg , AP gq , and 
APg g . They will be derived in the remainder of this section and agree with the limit x — > 
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of the leading order splitting functions of 0. In addition to these expressions for the 
kernels, we find the region of integration which gives the double logarithmic contributions. 
Next we will have to consider the bremsstrahlungs gluons. As the result of this section, we 
obtain all the ingreduents for the infrared evolution equations that will be discussed in the 
following section. Throughout this paper we will use the Feynman gauge and the Sudakov 
representation for the quark and gluon momenta: 

fci = antf + PiP + k Ti , (2.6) 

q' = q + xp, (2.7) 

and 

dki = -daidPidk^dipi. (2.8) 



For the beginning let us briefly review j|] the structure of the double log contribution 
in the non-singlet case. The imaginary part of the amplitude T 3 comes from the s-channel 
graph Fig.l, and the spin structure of (M L -M R ), corresponding to the quark line has the 
form 

tr [lsPl»{q + p)lv) = -lie^a^p 13 , (75 = -ieap^l^l 1 ! 5 ) ■ (2.9) 

which coincides with the tensor ie liva i3q a s 13 ■ m in eq.(2.2) for the longitudinal spin vector 
s /3 _ p/3 j m Q f foe initial quark. To obtain the leading logarithm in the first loop Fig.2a, 
one has to cancel one of the two transverse momenta factors k% from the t-channel quark 
propagators. In the trace 



tr 



IbPla'k-i^q + k)^ v k^ c 



2tr 



l^pkjnq'juk 



(2.10) 



we have to keep the largest momentum q' in the s-channel quark propagator (q + k), and 
retain only the transverse components of the photon polarization vectors (/i and v) (see 
eq. (2.9)). Therefore one can omit the aq^ component of k^ 



tr 



aq'q' 



aq tr [. 



as q 



rl 







(2.11) 



and write 



2tr 



l^pkl^qivk = 2k\tr [75^^ 



(2.12) 



Together with the quark propagators (1/kf) and the colour factor Cf = (N% — 1)/2N C it 
gives the logarithmic loop integral: 



2C F dkf , 2 
- ; — c^ s —— ■ k, 
4tt k 4 1 



(2.13) 
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The same reasoning is valid for the second loop fig. 2b ||. As a result, we obtain for the 
gluon rung inside the fermion ladder AP qq : 

AP qq = C F . (2.14) 

This result is also valid for the flavor singlet. 

Turning now to the singlet case, in the second loop approximation a new type of diagrams 
appears, contributions with t-cannel gluons. The first diagram of this type, Fig. 3, contains 
two fermion traces. The lower one, originating from the initial quark spinors, has the form 
eq.(2.9) with the gluon momentum (— k a ) instead of q a . It is antisymmetric in the polariza- 
tions of the two gluons with momentum k\. Thus we cannnot assign the same (longitudinal) 
polarizations to both gluons: otherwise we would have Mr = Ml and T 3 = \{Ml — M#) = 0. 
The best one can do in order to retain the largest power of 1/x (in the small x region) is 
to assign the longitudinal polarization oc q' , to one of the two gluon and the transverse 
polarization (e' v = e u r t ) to the other gluon. In the proton (quark) rest frame with k u = k x 
this means that e M / = e z , e u > = e y , and the lower trace 

tn = -tr (7 5 p<y£;i<v) = ^ie^^a^p 13 = -Aie z ^\p || k lt \e yv >. (2.15) 

To obtain the logarithmic integration over k± we have to find another factor ku from the 
upper trace, i.e. to take the component k\ t in the quark propagator (k 2 — ki) — the only 
one, which contains the momentum ku- The product of three 7-matrices e^, ki t e u i, in the 
upper trace may be written as 

e^'zkitxeu'y = -«757o|^it| • (2.16) 

Therefore the trace corresponding to the fermion loop in Fig.3 takes the same form (2.10) 
as for the Fig. 2b 



tr 



IsPhl^qivh , (2.17) 



and in this way the polarization of the target quark is transfered to the fermion loop. So the 
upper quark loop can be calculated in analogy with the one-loop graph in Fig.2a. As usual, 
the DL-contribution comes from the kinematical region of 

1 > A > (3 2 > . . . (2.18) 

and 

...<«!< a 2 <l. (2.19) 

where all the s-channel (horizontal lines) particles are on mass shell.0 In this case a^s <^ kf t , 
so kf w kf t . 



4 The ordering 1 > > f3i + \ > 0, < < a J+1 < 1 arises from the energy-momentum conservation, 
and within the DLA one can change it by the strong ordering eq.(2.13); this is the largest region which may 
give the DL-contribution. 
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Finally, the double logarithmic contribution of the diagram Fig. 3 to the g\ structure function 
is: 

*.~(?<M£)7£/¥t- 

From this expression we extract the results for the transitions AP gq and AP qg : 

AP gq = 2C F , (2.21) 
AP qg = -2T f . (2.22) 



It is interesting to note that the element corresponding to the gluon to quark transition 
(i.e. AP qg splitting kernel) — 2Tf is negative. From the physical point of view it means that 
the slow (low x) produced quark has the helicity opposite to the one of the parent gluon; 
the initial gluon helicity goes to the fast antiquark with larger x. 

At this point it may also be worthwhile to stress the change in the energy dependence of g\ 
compared to the unpolarized case, where the longitudinal polarization of both t-channel glu- 
ons leads to the behavior fi(x) oc 1/x, i.e. to a singularity at angular momentum j = 1. In 
the present case we cannot have two longitudinally polarized t-channel gluons and therefore 
loose one power in 1/x: 

gi (x) ~ x° ■ F(ln(l/x)) at x -> 0. (2.23) 

Consequently, the rightmost singularity of the function gi in the complex angular momentum 
j-plane lies at j = 0. The t-channel quark exchange gives the same power of 1/x at small x, 
and therefore quark and gluon structure functions mix at small x. This is in contrast to the 
unpolarized case, where at small x the gluons dominate by one power in 1/x. 

To complete our discussion of the rungs, we have to adress the diagrams with gluon rungs in 
a gluon ladder. Such a contribution appears first in the 3 d loop, and it is illustrated Fig.4. 
To simplify the kernel we will write the triple-gluon vertex as a sum of three pieces (Fig.5a): 

= g^u(h + k 2 ) p + g w {k 2 - 2k x ) v + g vp (ki - 2k 2 )^ (2.24) 

and note the following: 

1) due to the antisymmentry in the indices /iV the graphs 5b, c do not contribute to g\ 
(Fig.4); 

2) due to gauge invariance we have M^^ki^ = M^ u /ki u / = and N^i^ik^n = N^i^ik^n = 0; 

3) graph Fig.5d does not give the leading logarithm, as the product {k-y + k 2 ) 2 = 2{k\ + 
(note that the s-channel gluon k\ — k 2 is on shell and (ki — k 2 ) 2 = 0; the term ~ k\ (or k 2 ) 
cancels one of the gluon propagators and kills the logarithm ^ in the integral over k± t (or 
k 2t ) ). Thus the product of two triple gluon vertices in Fig.4 takes the form: 

H^'u'^'u" = 4( S^'n" k 2u /ki v " + Su'^'k^'ki^" — 5^ u »k 2u rki fM » — 5 u i^ki u ^k 2lJL ^ ) (2.25) 
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which corresponds to the diagrams Fig. 5c. 



Together with the colour coefficient N c it means that the insertion of an extra s-channel 
gluon into the gluon loop of Fig. 3 leads to the double log integration 



4Ar ^ up ^ (2 _ 26) 



a s r d(3 2 f dk 2 
; 2W J ~~kj 

Indeed, in order to keep the largest power of 1/x and to save the logarithm in the integrals 
over k t we have to choose in (2.25) the transverse components of k\ and k 2 and the longitu- 
dinal indices in 5 M y>, 6 v r v », 5 M v> or <^,y. Such a configuration conserves the main structure 
of the gluon loop fig. 3, inserting instead of the spin part of the transverse gluon propagator 
—e v ' y e v " y (in Fig. 3) the expression (e^k lt )(e L ,'k 2t ) = — \k u \\ k 2t \sm 2 Lp (here ip is the angle 
between the transverse momenta k\ and k 2 and we take into account the fact that e v > _L k u 
and e„» _L k 2t ). An additonal power of \ku\ and \k 2 t\ comes from the traces of the lower and 
upper quark loops (like in the case of eq.(2.15)) and (2.16)). So after the integration over 
the azimuthal angle ((sin 2 if) = 1/2) the vertex eq.(2.25) gives us the result (2.26) From this 
we extract the gluon rung AP gg : 

A gg = 4N C = AC A (2.27) 

Finally, let us collect our results for the four different rungs. We define a matrix M as 
illustrated in Fig. 6, which contains the splitting functions AP^-: 




(here Ca = N, Cp = 2 Z , and Tf — ^ are the usual SU(N) color factors; note that we 
have chosen to put the gluons into the first column and row) . This matrix will be used in the 
following section where we shall derive the infrared evolution equations. For later purposes 
it will be convenient to consider also the color octet t-channel. In this case the color matrix 
analogous to (2.28) reads: 



From our previous discussion we have also obtained the general pattern of the region of 
phase space which gives the double logarithmic contributions. The limits of integrations fol- 
low from the ordering condition given in eqs.(2.18),(2.19). In terms of $ and k ti it means 
that 



15 



0i 



(here k 2 ti > 0). 



(2.30) 
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k 2 ■ 

Indeed, let us first take fa + \ <C fa and kt j+i > fcti. In this case a.i+\ ~ t '! j +1 ^> a 



we have used the fact that s-channel partons with the momenta hi — ki + i are on mass shell; 

1,2 +1 — ai + ifas = 0). On the other hand, if kt,%+i < k t ,i one has a^+i ~ j^. 

t+1 ~ "'t.i+l' 



(hi — fcj+i) ~ — k 2 i+l — ati+ifas = 0). On the other hand, if kt,%+i < k t)i one has a^x 
vc 

a i+ iP i+ iS <C kf i+1 . In other words, our condition looks as 



In order to save the leading logarithm, we have to satisfy the condition: k 2 +l rs —k 2 i+1 , i. e. 



k 2 Ifi 

fa+l s ftiS 



a l+1 = -£±i » = a, (2.31) 



and as in P we can simply use the ordering eq.(2.19) with the oti (eq.(2.31)) instead of the a«. 

In the final part of this section we have to consider nonladder diagrams as illustrated in 
Fig. 7. We will call a non-ladder gluon "soft" if its transverse momentum is smaller than 
the momenta of all the partons comprised by the non-ladder gluon. According to Gribov's 
fcr-factorization theorem 0] [] the whole amplitude of the 'soft' gluon emission can be writ- 
ten as the bremsstrahlung from one of the "external" lines of the block comprized by this 
gluon. Therefore the double logarithmic contribution coming from such a 'soft' gluons can 
be summed up with the help of the infrared evolution equation |§, in the same way as it was 
done for the non-singlet structure function gi in[|J. When summing over all possibilities of 
attaching the soft gluon to the external legs, we get a total color factor which depends on 
both the total t-channel color quantum number and the type of incoming partons. We will 
need the color singlet channel. For incoming gluons and fermions the color factors are Ca 
and Cf-i resp. In matrix notation we define 



Go = / n • (2-32) 




Finally we note that non-ladder gluons with k' t larger than the transverse momenta in the 
part of the ladder, which is comprised by them (they will be called 'hard') do not give double 
logarithms. For a nonladder gluon that runs across the ladder from one side to the other 
(e.g. from the lower left to the upper right leg), the large momentum has to flow through 
some internal small k ti ladder propagators, and the large momentum k' t changes the normal 
factor to l/k' 2 , in this way killing the leading logarithm dk\lk\ (see Fig. 7). These hard 
nonladder gluons therefore do not contribute to the double logarithmic approximation. Next 
we consider vertex correction (Fig. 7b). Then we can say that in the ultraviolet (large k' t ) 
region there are no any double logs in the Feynmann gauge for the vertex function. So, one 
can anticipate that there are no DL-correction coming from large- k t - region in the Feynmann 
gauge for the vertex function. However, we have to be more carefull here. For the unpolar- 
ized case there exists an example (BFKL) where a special cutting of the vertex-type diagram 
does give rise to a double log contribution. Indeed, let us consider the non-ladder on-shell 
gluon k' added to the amplitude: the loop integration yields , two propagators - 

For QCD the Gribov's theorem was considered in more detail in §, § and |@ 
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1/q 2 ~ ^7^; 1/^1 ~ V^it an d the spin part of the propagators, which for the nonsense, 
longitudinally polarization (g'" u , « Q' u Pv> /(Q'p)) in the DL kinematical region /3 2 -C /?' — /3i 
(a g ^> a') gives the factor ~ — (qp)(2kiQ')/(pQ) ~ a q (3's. This factor cancels the propagator 
1/q 2 , while the next propagator 1/A;f t ~ provides the logarithmic integration dk l2 /k' 2 in 
the region fc£ > k u . As it is known ]TJ], the sum of such contributions is equal to the ladder 
contribution (Fig. 4) but has an opposite sign. So it cancels the double logs coming from 
k' t > k ti and restores the conventional DGLAP ordering k tt i + i > k t>i for the gluon loops. 
Fortunately, this does not happen for the spin dependent structure function g-±. In this case 
the vertex of the "nonsense" gluon k' emission changes the sign, if unstead of the longitudi- 
nal polarizations of i-channel gluons k\ and ki (in the left side of Figs. 7c, d we consider the 
transverse polarization e\ _L k u ^\. Thus the amplitudes of the Fig.7 type, where the gluon 
k' is emitted one time by the transverse gluon and another time by the longitudinal one, 
cancel each other and we come back to the ladder configuration with the ordering eqs.(2.30), 
(2.31). 



3 Infrared Evolution Equations 

In this section we construct the infrared evolution equations which are necessary for the 
calculation of T3 and g±. We shall follow [|J, and we begin with the amplitude T3 which, 
following the discussion of the previous section, consists of the two components: 

(from now on it will be understood that we consider the singlet part only, and we suppress 
the subscript "S".) The structure function then follows from the relation (2.3), and we 
have to take into account both DL-contributions and in -terms. We write T3 as a Mellin 
transform: 



too 



du 



T 3 = I . —(^r^)R(u,y), (3.2) 



voo 



where R(u>,y) is a two-component vector, defined on analogy to (3.1). The signature factor 
f M = is: 

= €ZZ1 „ (3.3) 



6 To save logarithm in d/3' / j3' integration the gluon k' should be the longitudinal, nonsence one; on the 
other hand the polarization vectors of the ^-channel gluons ki (or £2) and k\ should be — one transverse 
and one longitudinal (for the small-x limit of g±) as it was discussed just before eq. (11). If the vector e\ || p 
than the leading contribution comes from the graphs Fig.7e,f and T^ up ■ p^Q'^Qp/ (pQ 1 ) — ~^jr = — a ii 
while for the transverse vector e\ — en — en only the diagram Fig.7e do work and r \ V pZ\tii&itvQ ' p / 'pQ' ~ 
+ai + 0(a 2 ); (a 2 -C Qi). 
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and 

y = H%). (3.4) 

Thus 

n OR , d , _ . . 



Eq.( p.5| ) represents the left-hand side of the IREE for T 3 which is illustrated in Fig. 8. 
The right hand side is obtained from the observation that the dependence upon the cutoff 
H resides in the intermediate state with lowest virtuality (Fig.8): the /^-derivative of the 
amplitudes are equal to R times quark or gluon scattering amplitudes with the external legs 
having transverse momenta close to /i. We are thus lead to the definition of a quark-quark 
scattering amplitude for which we use an integral representation of the form (3.2). The 
partial wave will be denoted by F qq . More general, we introduce the four amplitudes F qq , 
F qg , F gq , and F gg , and in analogy with (2.28) we combine them into the two by two matrix 
F 

Fo =(f 9 f") (3 - 6) 

(it is the analogue to /q m 0)- 111 terms of this Fq, the vector evolution equation for R 
becomes (Fig.8): 

In analogy with [|| || the evolution equation of Fq has the form (Fig. 9a): 

Q 2 Q 2 1 

uj z7t z uj z Sn^uj 

Here we have used the the matrices M and Go defined in the previous section. The second 
term on the rhs of (3.8) corres ponds to the gluon bremsstrahlung diagrams: in analogy to 
the matrix Fq which carries color zero we define the matrix F$ of color octet amplitudes (it 
is the analogue to in [[|). These amplitudes satisfy the evolution equation similar to 
(3.8): 

F 8 = 9 -M s + P^^Lf 8 (cu) + ^F 8 (uj) 2 . (3.9) 
uj Stt^uj duj &1T Z UJ 

The matrix M 8 is taken from the previous section, and the color factor Ca in front of the 
second term on the rhs is the analogue of the matrix Go in (3.8). The difference between 
Ca in (3.9) and G in (3.8) is due to the fact that for the positive signature amplitude F 8 
the sum of the two bremsstrahlungs diagrams (illustrated in Fig. 9b) is independent of the 
type of the incoming partons, and the matrix of color factors G$ becomes Ca times the unit 
matrix. 



F (u) = y -M - -JL-^GMu) + —F {uj) 2 . (3i 
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4 Solution of the evolution equation 

The solution for g\ is obtained by solving eq.(3.9) for F&, then eq.(3.8) for F and eq.(3.7) 
for R (the latter makes use of the Born approximation Rb as initial condition for R), and 
finally inserting R into (3.2). The final answer for T 3 is the two component vector: 



Let us go through these steps in somewhat more detail. 

We begin with the equation for F%. We first diagonalize the Born term, i.e. the matrix 
Mg. The eigenvalues are: 

^ = ' 2 1/2iV ± \4Wa - 1/2NY - AC A T f . (4.2) 

Let and denote the two eigenvectors of M 8 and E 8 = (eW,e^) the matrix composed 
of them. Then 

M 8 = E 8 M 8 E^, (4.3) 

where M 8 = diag(\ 8 , A 8 ). Consequently, eq.(3.9) becomes diagonal if we transform to F 8 : 

F 8 = E 8 F 8 E^. (4.4) 
The solutions for the two components of F 8 are: 

Fs* = Ng 2 -?-(\ne* 2 /*D p± (z)), (4.5) 
where D p denotes the parabolic cylinder function with 

A (±) 



and 



z = —, luo = y/NgySn*. (4.7) 



With this solution for F 8 we return to the evolution equation (3.8) for F which is solved by 
the (matrix-valued) expression: 



4^ F ° ~ 1 " \ 



2(nu) 



;M + 



9 Z 



tG F 8 



(4.8) 
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This has to be substituted into the final formula (3.2) for T 3 . The leading behavior at small 
x comes from the right-most singularity in the uo plane. Similar to the nonsinglet case, this 
singularity is due to the vanishing of the square root in (4.8), i.e. we need to determine 
the zeroes of the eigenvalues of the matrix under the square root. Similar to the discussion 
after (4.2), the diagonalization of this matrix is done by the matrix E consisting of the two 
eigenvectors and of the matrix under the square root. 

Let us first return to the eigenvalues of M 8 determining the asymptotic behavior of F s . 
For four active flavors the eigenvalues are: 



\f = 2.92 ±1.88= -° . (4.9) 




Neglecting fermions one obtains A = 6, whereas the pure fermionic case corresponds to 
A = — 1/6. Obviously, the larger of the two eigenvalues is not too far from the gluonic case. 

The accurate values have to be found from a numerical computation of the parabolic cylinder 
function. Before we quote the results of such a calculation it may be instructive to discuss a 
few approximations. Turning to the square root expression in (4.8), we are searching values 
of uj (or z) for which one of the two eigenvalues of the matrix under the square root in (4.8) 
becomes zero. To obtain a first estimate, we simply neglect the term containing F s (the 
exact calculation will show that this approximation is not too bad): then we can diagonalize 
F by simply diagonalizing M . The eigenvalues are 

A± = \{AC A + C F ) ± \4Wa + Cf) 2 - 16C A C F - 16C F T f (4.10) 
with the corresponding eigenvectors: 

« (+, = (j«). * (+, = ^ = ^ ( 4 - u ) 

For the diagonalization we need the matrix E = (e^ + \ e^) and its inverse 

E, 1 = - \ ( \ "f V (4.13) 
u 1 - x+x- \ -x + I J 

The two values of z are (for four active flavors): 

(4.14) 

(the pure gluonic and fermionic cases give z = 4 and z = 4/3 resp.). With the eigenvalues 

(4.15) 
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the matrix E is found to be: 



For the partial wave R(u,y) in (3.2) we obtain: 

Retaining in F only the leading upper component, we find for the behavior of R near the 
square root branch point at cu = u s : 

The fact that the matrix elements in the rightmost column are both negative has the impor- 
tant consequence that the leading contribution at small x changes the sign relative to the 
input distribution. Finally, 



8v^F (ln(l/x))3/2 (A ^' AS )^> « )(^) ^ + ^HuT/^ 



where the vector R is from (4.18), and the (transposed) vector (Ag, AS) represents the 
initial conditions of the gluon and quark polarized distributions. It should be kept in mind 
that here we have retained only the leading singularity uo s of (4.14). With ujq from (4.7) and 
a s = 0.18 we find uj s = z s ^a s N c /2n = 1.12. 

Finally, taking into account also F$ we quote the result of a numerical calculation. For 
the larger of the two z- values we find (n/ = 4): 

z s = 3.45 (4.20) 

(the pure gluonic case would have given z = 3.66). With c^o from (4.7) and a s = 0.18 we 
find 



,yJa,N c /2Tr = 1.01, (4.21) 
and with x + = 0.29, x~ = 0.43 

» f . ( hU -°- 50 \ 2 (Q 2 Y s/2 ( o \ 

-0.14 J ^ J U^J' ( } 
i.e. the negative sign persists. The result for g 1 becomes: 

M = ^4±i^! (A9 , AE)Ji(w ,,, )( I r ( 1 + g™)) (4-23) 



8^ (ln(l/x)) 3 / 2 v y ' ' v s,y; V V v lnl/x 
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with R from (4.22). 



Let us finally compare our result with the fixed order calculations. The region where we 
expect our result to coincide with the fixed order calculation is ^fo~ s <w<l. We therefore 
expand our anomalous dimension matrix in powers of g 2 /uj 2 . To leading order we obtain: 

It agrees (as we have already mentioned before) with the singular part of 0. The 
next-to-leading-order matrix has the form: 

m _ a 2 1 / 32C1 - l6C F T f -16C A T f - 8C F T f \ 
75 ~ [ 4^ } ~z{ 16C A C F + 8C F AC 2 - \QC F T f + ) ' ^ 25j 

After transformation to xb, it agrees with the leading powers on In x given in eq.(3.65)-(3.68) 

of nrai. 



5 Discussion 



The main result of this paper is the power-like behavior of the flavor singlet part of g\ at small 
x (eq.(4.23)). The (leading) power is by a factor 2.6 larger than in the nonsinglet case. This 
effect is mainly due to the t-channel gluons states, which have a much larger color charge 
(cf.eq.(2.28), (2.29)) than the quarks (comparing the pure gluon and the pure fermionic case 
and neglecting the nonladder gluons one finds u) smglet / ' 0J nonsin 9^t _ ^J4C A /C F = 3). As it can 
be seen from a comparison of (4.14) and (4.20), the influence of the bremsstrahlung gluon is 
of the order of 10%. Comparing the flavor singlet with the nonsinglet, one notices that the 
nonladder bremsstrahlungs gluons act in quite different ways: in the present case, due to 
the positive sign of the largest matrixelemnt of M 8 (M 8 n ), the leading cu-plane singularity 
moves to the left, i.e. the exponent of 1/x decreases from uj s = 1.12 to uo s = 1.01. For the 
flavor nonsinglet, on the other hand, the nonladder gluons lead to an increase of the exponent. 



As to the polarization of the t-channel ladder gluons, the following pattern has emerged. 
In the unpolarized case, both t-channel gluons are longitudinally polarized (i.e. they are 
in the nonsense helicity state), and the resulting behaviour is g(x,Q 2 ) ~ (1/ 'x) l+0<yaa \ In 
the polarized case, this leading contribution of gluon polarizations cancels, and the next- 
to-leading configuration appears: one gluon is still longitudinally polarized, but the other 
one has transverse polarization. The small-x behaviour is now of the form gi(x.Q 2 ) ~ 
(1/ 'x) 0<y ^^\ If both gluons are transversely polarized, the small-x behaviour is further sup- 
pressed: ~ (l/x)~ 1+0( -^*°\ This contribution has been studied in [O . 



Compared to the small-x prediction of the standard GLAP evolution equation, the main 
reason for the strong enhancement of our calculation lies in the different regions of phase 
space. Whereas in the GLAP case we have strong ordering in transverse momenta of the 
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ladder partons, our double logarithmic approximation has a much larger region of integra- 
tion. A good way of illustrating the difference is the diffusion in In k\ which first has been 
observed for the BFKL |T3j Pomeron. When going, within ladder graphs, from one rung "i 



to the next rung "i+1", ordering in transverse momenta means kxi < kn+i- In the BFKL 
approximation, the distribution in In kn+i is a random walk distribution, i.e. when climbing 
up the ladder up to rung "i+1 ", we either add or subtract, at each step, a fixed steplength A. 
In the present double logarithmic approximation, on the other hand, the steplength is not 
fixed but can vary between zero and lnl/x. From this qualitative argument it is clear that, 
in the polarized case, the difference between the GLAP small-x behaviour and the double 
logarithmic prediction should be larger than, in the unpolarized case, the difference between 
GLAP and the BFKL behaviour. 

Finally, at least a few words should be said about the phenomenological applications of 
our calculation JT5], [H| [TTJ. Of particular interest is, of course, the question at what values 



of x and Q 2 the GLAP extrapolation should be replaced by the result of this paper. Leav- 
ing a more detailed discussion for future publications, we only want to make a few general 
remarks. In this paper we have derived expressions for Greens functions - G qt9 (x, x'; Q 2 , fi 2 ) 
- which describe the evolution from the initial point (a/, fj, 2 ) to the final point (x, Q 2 ). If the 
initial polarized quark and gluon distribution do not exhibit any particularly strong variation 
in x, we can start from the approximation: 

gi(x, Q 2 ) = AE(x, /i 2 )G g (x, x; Q 2 , /i 2 ) + Ag(x, fi 2 )G g (x, x; Q 2 , fi 2 ) (5.1) 

where x denotes the mean x-value of the initial distributions. As a general feature, our 
Greens functions rather strongly depends upon the momentum scale fi 2 . 

Our analytic expression (4.23) represents only the leading term of the full amplitude (4.1). 
We have checked numerically that, as a function of x, the full amplitude (4.1) quickly shows 
the power behaviour predicted by the leading approximation (due to the factor lnl/x there 
is, however, some difference in the overall normalization). Furthermore, as (4.22) suggests, at 
small x the gluon will dominate over the quark component. A numerical calculation shows 
that this, again, happens rather rapidly: the ratio of gluon and quark Greens functions 
Gg/Gg approaches its asymptotic value 3.45 approximately at x — 10~ 3 . 

An important aspect of our result (4.23) is the sign structure. Our formula predicts that, 
at small x, the sign of g\ should be opposite to that of the initial distributions. The ex- 
perimental data Q indicate that at small x the proton function g\[x) is positive, and even 
the neutron function g\ for x < 0.01 shows the tendency to change the sign and become 
positive. So, at first sight, there seems to be a problem. A possible way of resolving this 
conflict could be the following. Whereas the value of AS is known (and positive), we do not 
know the initial value of Ag. In ref.|]18[ it was shown that in a resonable model of the nu- 
cleon (non-relativistc quark or bag-model) this value may be negative, of about Ag ~ —0.4. 
Taking into account that at small x the gluon Greens function will dominate over the quark 
component, the structure function gi may very well become positive - both for the proton 
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and the neutron - and thus be consistent with the behavior observed in the data. 



All these arguments are clearly very qualitative, and a more detailed quantitative analy- 
sis is of great importance. For the nonsinglet case, a first attempt has been made in |15| 
to estimate numerically how large the deviations of our double log formula from the usual 
GLAP predictions might be. Clearly any such estimate will strongly depend upon the input 
distribution. In [15| it was concluded that the deviations may be very large (up to an order 
of magnitude) but the difference between the nonsinglet gi(x,Q 2 ) and fi(x,Q 2 ) will not 
exceed 15 % within the HERA kinematical range. A somewhat different analysis has been 



performed in [17]]. The situation with the singlet contribution will be more involved: not 
only is the power generated by the double logs larger than in the nonsinglet case, also the 
sign structure of the leading term (as mentioned in the previous section) will make the result 
more dependent upon the input. Clearly a careful study in this direction is strongly needed. 
On the longer run, also an improvement in accuracy of our double logarithmic approximation 
should be persued. 
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Figure Captions: 

Fig.l: Born approximation for M s . 

Fig. 2: First corrections to M . gluon rungs in the fermion ladder. 

Fig. 3: Lowest order contribution with t-channel gluons. 

Fig. 4: The first gluon rung. 

Fig. 5: Calculation of the gluon rung in Fig. 4. 

Fig. 6: Matrix notation in Eq.(2.17). 

Fig. 7: Nonladder graphs. 

Fig. 8: Illustration of the infrared evolution equation (3.7). 

Fig. 9: (a) Illustration of the evolution equation (3.9); (b) origin of the color factor in front 
of the second term in (3.9). 
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